Dark Matter Spikes and Annihilation Radiation from the Galactic Center 
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The annihilation rate of weakly interacting cold dark matter particles at the galactic center could 
be greatly enhanced by the growth of a density spike around the central supermassive black hole 
(SBH). Elere we discuss the effects of hierarchical mergers on the central spike. Mergers between 
halos containing SBHs lead to the formation of SBH binaries which transfer energy to the dark matter 
particles, lowering their density. The predicted flux of annihiliation photons from the galactic center 
is several orders of magnitude smaller than in models that ignore the effects of SBHs and mergers. 
Measurement of the annihilation radiation could in principle be used to constrain the merger history 
of the galaxy. 



The dark matter puzzle is one of the central challenges 
facing particle physics and cosmology . A popular can- 
didate for non-baryonic cold dark matter (CDM) is the 
lightest supersymmetric particle, plausibly the neutralino 
X ^ . The mass of the neutralino is constrained by ac- 
celerator searches and theoretical considerations of ther- 
mal freeze-out to lie in the range 30GeV < < lOTeV 
[Q, ||. Neutralinos are generically found to decouple at 
a temperature that is roughly M^/20, which means that 
they are nonrelativistic already at decoupling and hence 
behave like CDM. 

Dark matter particles may be detected directly, via 
laboratory experiments or indirectly, via their an- 
nihilation products 0. Indirect schemes are typically 
based on searches for gamma rays from neutralino an- 
nihilations in the dark matter halo of the Milky Way 
(MW) galaxy ||, |, Since the photon flux depends 
on the squared density of neutralinos integrated along 
the line of sight, the signal is greatly enhanced in direc- 
tions where the dark matter is clumped. This includes 
the galactic center, where the density in a smooth halo 
would be maximum, as well as any lines of sight inter- 
secting the centers of relic halos that orbit as subclumps 
in the MW halo The signal from the galactic 

center is further enhanced if there is a CDM "spike" as- 
sociated with the central supermassive black hole (SBH). 
Adiabatic growth of a SBH at the center of a pre-existing 
halo produces a power-law distribution of matter around 
the SBH, a "spike," with density p - r"'', 2.2 < 7 < 2.5, 
r < rh = GM,/a'^ M, is the SBH mass and a is 
the ID velocity dispersion of the dark matter particles. 
For the MW SBH, M. ~ 2 - 5 x IQ^Mq and 
r/j ~ 1 pc. These spikes lead to predictions of higher- 



than-observed rates of 7-ray annihilation products! 17 



One element missing from earlier discussions of dark 
matter spikes is the destructive effect of hierarchical 
mergers. A dark matter halo as massive as that of the 
MW, M ~ IO^^Mq, has almost certainly experienced a 
significant merger event since a redshift of 2; ~ 2. Fur- 
thermore, SBHs with masses of ~ IO^Mq were present in 
at least some halos already at redshifts of 5 — 6 18, p^, 
and SBHs probably acquired most of their mass by a red- 



shift of 2 — 3, the epoch of peak quasar activity P0|. In 
the CDM paradigm, big halos grow through the buildup 
of smaller ones if more than one of the progenitor 
halos carried a central SBH, a binary SBH will form fol- 
lowing the merger . Formation and decay of the SBH 
binary transfers energy to the background particles and 
lowers the density of matter in a region within a few x rh 
around the binary 1 23 , 2^ , roughly the scale of the spike, 
before the SBHs coalesce due to emission of gravitational 
radiation psf . In this paper we compute the dark matter 
density profiles of merging CDM halos containing SBHs 
and show that the net result of including the SBHs is to 
substantially lower the predicted flux of annhilation pho- 
tons compared with halos lacking SBHs. Our results are 
relevant also to models that relate the origin of ultra high 
energy cosmic rays to the decay of unstable superheavy 
relic particles in the galactic halo |2^ . 

We first calculated the probability that the MW halo 
has experienced a major merger event by generating 
multiple realizations of the merger history of a halo of 
lO^^AfQ using the algorithm described by Somerville and 
Kolatt , as applied to the standard ACDM cosmolog- 
ical model (f2o = 0.3, Aq = 0.7). This algorithm accu- 
rately reproduces the merger histories of halos seen in A^- 
body simulations of structure formation |^ although at 
a value of the scale factor that may differ as much as 20% 
from the true value |29 . We recorded all merger events 



such that the mass of the smaller of the subclumps in- 
volved in the merger was above a given limit M\i^. From 
600 of these realizations, we then computed the proba- 
bility that a halo of IO^^Mq at redshift z — Q had in its 
merger history at least one such event since a redshift of 
z = 2. Fig. 1 shows this probability as a function of 
Afiim- With 68% confidence, a halo the size of the MW 
has experienced at least one merger since z — 2 m. which 
the mass of the smaller of the merging halos was above 
2 X 1O^^M0, implying a mass ratio between subclumps 
of 4 : 1 or less. Mergers with mass ratios of 9 : 1 or less 
occurred with greater than 90% probability. 

To investigate the detailed effects of mergers on dark 
matter spikes, we used an A'^-body code to simulate in- 
teractions between CDM halos containing central SBHs. 
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TABLE I: Parameters of the iV-body integrations. Mi (M2) 
is the mass of the large (small) halo; p is the central halo 
density before growth of the SBH; M. is the mass of the SBH 
particle; N is the combined number of DM particles in both 
halos; aflnai is the final separation of the BH binary. 



Run 


M1/M2 


Pl/P2 


M./M 


N 


aflnal(pc) 


A 


1 


1 


0.01 


10^ 


0.91 


B 


1 


1 


0.03 


2 X 10^ 


0.71 


C 


3 


1 


0.03 


4 X 10* 


0.31 


D 


3 


1/3 


0.03 


4 X 10* 


0.42 


E 


5 


1 


0.03 


6 X 10* 


0.36 


F 


10 


1 


0.03 


1.1 X 10^ 


0.90 



iV-body models were generated from the spherical den- 
sity law 



PBuix) = Pox ^-^(l + x^-^) \ 



r/ruM, (1) 



the "Moore profile" Q i '^^^ ^ set of profiles found to 
accurately describe the halos generated in CDM struc- 
ture formation simulations; r^u defines the radius of 
transition between the inner cusp and the steeper outer 
falloff. Our choice of the Moore profile is conservative in 
the sense that other proposed fitting functions |3|] 
have shallower central density cusps implying less promi- 
nent spikes. The Moore profile has a divergent mass; 
we truncated it spatially and generated velocities of the 
dark matter particles from a distribution function 
which generates an isotropic velocity distribution near 
the center and increasingly circular orbits approaching 
the truncation radius. This procedure guarantees a state 
of detailed dynamical equilibrium in spite of the model's 
hard edge. Since the dependence of the central densi- 
ties of CDM halos on halo mass is not well determined, 
we investigated two values of the central density ratio 
Pi/P2 = {1,1/3} for one of the unequal-mass mergers, 
where p is defined as the asymptotic central density of 
the halos before growth of the SBH. A "black hole" was 
inserted into each iV-body halo by slowly increasing the 
mass of a (Newtonian) central point particle from to 
M, (see Ref. [|4|). We verified that the resulting den- 
sity profile satisfied pdm ^ r~^-^, r < r/j as predicted 
by adiabatic theory ([|l^; Fig. 2). Since the scale length 
of the MW halo is believed to be tdm ^ 20 kpc 
the dark matter density profile is essentially scale-free, 
Pdm ~ r~^'^, at all radii of interest and the mass chosen 
for the particle representing the MW SBH is fairly arbi- 
trary. We tried two values, M,/Mdm = {0.01,0.03}; the 
larger value was preferred since it resulted in a lower am- 
plitude of Brownian motion of the SBH particle . The 
same mass ratio was adopted for the SBH in the smaller 
halo. The parameters of the A'^-body runs are given in 
Table |. 

The mass scaling of our models was fixed by the mea- 
sured mass of the MW SBH 111, fsl. RE; we assumed 
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FIG. 1: Probability that the merger tree of a CDM halo of 
mass 10^^ M0 at 2: = contains at least one merger event 
between z — and 2 = 2 in which the mass M2 of the smaller 
subclump exceeded Mum. 



M, = 3 X IO^Mq. The length scale was determined 
by assuming that the dark matter in the MW produces 
a circular rotation velocity of 90 km s^^ at the solar 
circle, Rq — 8 kpc and that its density profile is 
PDM c>c r~^-^ inward of Rq. 

Simulated mergers of the halo models were carried 
out using a new, general-purpose iV-body algorithm [3^ 
that combines the elements of a hashed tree code with 
a quadrupolar expansion of force moments for the bulk 
dynamics, and the highly accurate Hermite predictor- 
corrector scheme for near-neighbor and massive-particle 
interactions. The code implements individual block time 
steps, individual stepping and softening criteria, and full 
functional parallelization. Interactions between SBH and 
dark-matter particles were unsoftened. Calculations were 
carried out using 16 processors on the Rutgers Sun HPC- 
10000 computer. Integrations were terminated when the 
separation between the SBH particles was less than 1 pc 
(see Table 1). 

Fig. 2 shows the effect of mergers on the dark-matter 
density profiles. The steep, p ~ r"^ "' density spikes are 
destroyed in each case by the transfer of kinetic energy 
from the binary SBH to the dark matter particles. (Merg- 
ers without central SBHs tend to preserve spike slopes 
[24, |39|.) The energy transfer takes place in two stages 
[ 24 1 : before the two SBHs form a binary system, dynam- 
ical friction acting on the individual SBHs causes their 
orbits to decay; and after formation of a bound pair, dark 
matter particles which pass the binary within a few times 
the binary's semi-major axis are ejected by the gravita- 
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FIG. 2: Density profiles of tiie merged dark matter halos. 
Tiie origin is defined as tlie center of mass of tlie binary SBH; 
spfierical symmetry is assumed. Tfie curves labelled M and 
S are the density profiles of the large halo before and after 
growth of the SBH. 
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FIG. 3: Dimensionless integrated flux of the merger models 
as a function of the angular acceptance of the detector Ail 
(Eq. 4). The curve labelled M is the flux predicted by a 
Moore-profile halo without a SBH- induced spike (Fig. 2). 



tional slingshot. The result is a lowered density out to a 
"core" radius Tc of ~ 10 — 100 pc and a density profile 
that rises inward of Vc as a weak power-law, p ^ r~^'^ . 
The amount of damage done to the pre-existing spike in- 
creases with the mass of the secondary SBH, consistent 
with the expectation that the mass ejected by the binary 
is of order a few times M2 ||4^. Final mean dark-matter 
densities within 100 pc are 5 — IOM0 pc"'^ for all of the 
runs, rising inward to p 10^ Af© pc^'^ for the 1 : 1 
mergers and ~ lO^M© pc^"^ for the 10 : 1 merger. 

The differential photon flux along a direction that 
makes an angle "0 with respect to the galactic center is 

where fl is the solid angle, p is the neutralino density, and 
av is the annihilation cross section (independent of v for 
nonrelativistic particles) ; the sum is over all annihilation 
channels. = 2 for xx 77 and — 1 for xx ~* ^1 
Ip^ . We are principally concerned with the final line- 
of-sight integral, J{ip) — J^p'^dl, which contains all of 
the information about the halo density profile. Following 
earlier authors [Q, |l^, we write J in dimensionless form 
as 

= TTT^ ( / P^dl (3) 

8.5kpc \^0.3GeV/cm7 

where the normalizing factors for length and density are 
roughly the radius of the solar circle and the local density 
of dark matter respectively. Finally we average the flux 



over the field of view assuming a circular aperture of size 
Ari centered at -0 = 0: 

Fig. 3 plots (J) as a function of Afi at the end of the A^- 
body simulations. For a typical atmospheric Cherenkov 
telescope angular acceptance of Afl = 10~^ sr, (J) ranges 
from ^ 10^'^ for the 1 : 1 merger, to ~ 10^-^ for the 
3 : 1 mergers, compared with ~ lO'''^ for the initial halo 
model without a SBH-induced spike (computed assuming 
an inner cutoff at the Schwarzschild radius of the SBH). 
Thus the addition of SBHs to CDM halos results in a net 
decrease in the annihilation flux compared with SBH-free 
models, if mergers are taken into account. The values 
of (J) in Fig. 3 are nevertheless large enough to allow 
testing of large parts of the neutralino parameter space 
using instruments like GLAST Q. 

We note that the predicted flux is a strong function of 
the merger parameters when observed with Ail < 10"^ 
(Fig. 3), corresponding roughly to the sphere of influ- 
ence of the MW SBH. This fact might allow the merger 
history of the MW to be inferred from measurements of 
the annihilation flux on different angular scales. 

The addition of SBHs to CDM halos could result in 
even lower densities than shown in Fig. 2. Halos as 
massive as that of the MW are believed to have formed 
through a succession of mergers, and the damage done by 
binary SBHs would be to some extent cumulative, result- 
ing in shallower central proflles than found here p2| . 
In addition, binary SBHs may eventually coalesce due to 
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emission of gravitational radiation |2^]. The radiation 
so emitted carries linear momentum leading to a recoil of 
the SBH at a velocity Vr^con 10^ - lO^ km/s || ||, || , 
and possibly even higher if the SBHs were rapidly spin- 
ning prior to coalescence |46j . Recoil velocities of this or- 
der would eject the coalesced SBH from the nucleus and 
its subsequent infall would displace dark matter particles 
(e.g. Q). Quantitative evaluation of this effect will 
require more accurate estimates of Wrocoii based on fully 
general-relativistic calculations of black hole mergers. 
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